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ABSTRACT 

We show that the shear modes in the neutron star crust are quite sensitive to the 
existence of nonuniform nuclear structures, so-called "pasta" . Due to the existence 
of pasta phase, the frequencies of shear modes are reduced, where the dependence of 
fundamental frequency is different from that of overtones. Since the torsional shear 
frequencies depend strongly on the structure of pasta phase, through the observations 
of stellar oscillations, one can probe the pasta structure in the crust, although that is 
quite difficult via the other observations. Additionally, considering the effect of pasta 
phase, we show the possibility to explain the all observed frequencies in the SGR 
1806-20 with using only crust torsional shear modes. 
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The soft gamma repeaters (SGRs) are considered 
as one of the most promising candidate of magnetars, 
which are neutron stars with strong magnetic fields 
( Duncan fc Thompsonlll992l ). The sporadic X- and gamma- 
ray bursts are radiated from SGRs, while SGRs rarely emit 
much stronger gamma-rays called "giant flares" . Up to now, 
at least three giant flares have been detected, which are the 
SGR 0526-66 in 1979, the SGR 1900+14 in 1998, and the 
SGR 1806-20 in 2004. For each giant flare phenomenon, one 
can observe a decaying softer part (tail) for hundreds of sec- 
onds after the initial short peak in the hard part of the spec- 
trum. Through the timing analysis of these decaying tail, the 
quasi-periodic oscillations (QPOs) have discov ered, which 
are in the range from tens Hz up to a few kHz (|lsrael et al.l 
l2005l : I Watts fc Strohmaved 12009 ). Since the QPOs are be- 
lieved as the outcomes of the neutron star oscillations, the 
observations of QPOs in SGRs could be first evidences to 
detect the neutron star oscillations directly. 

The current understanding about the observed QPO 
frequencies in SGRs is as follows; a part of lower frequen- 
cies such as 18, 26, and 30 Hz in the SGR 1806-20 are as- 
sociated with the magnetized fluid core, while the others 
are interpreted as torsional shear modes of the solid crust. 
However, it seems to be more complicated to understand 
the oscillations of magnetized neutron stars completely. 
The oscillat i ons in the core become t he Alfven continuum 
(lLevinll2006l.l2007l:ISotani et al.ll2008al : IColaiuda etai]|2009t 
ICerda-Duran et al l 12009^ 7 and one might need to consider 
the coupling between oscillations in the fluid core and in 
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solid crust if the m agnetic field is stronger jLevinl |200d . 
120071 ). Additionally, Ivan Hoven fc Levin! (|201ll ) described 
how crustal modes may survive in the gapes left in the 
Alfven continuum. In the spite of these complexities, the 
crustal torsional modes can still emerge globally if the mag- 
netic field is not so strong, e.g., B < 10 14 G, and the fre- 
quencies of torsional shear modes are almost same values 
as those in the case without mag netic filed (|Gabler et al.l 
l201ll : IColaiuda fc Kokkotasl l201ll ). Thus, comparing the 
analysis of shear modes with the observed QPO frequen- 
cies, one can know the proper properties of the crust in 
neutron star. In fact, via the observed QPO frequencies, 
it could be possible to constra in on the stellar properties 
l|Samuelsson fc Ande rsson 200J ) and th e distribution of stel- 
lar magnetic field ( Sotani et al l2008bl ). 



It is considered that the neutron star crust exists from 
the bottom of the ocean of melted iron at a density ~ 
10 6 — 10 8 g/cm 3 inward to the boundary with the inner fluid 
core at a density of order the saturation density of nuclear 
matter p s ~ 3 x 10 14 g/cm 3 . Although nuclei in the crust 
form a bec lattice due to Coulomb interactions, according 
to the recent studies, the nuclear structure in the bottom of 
crust could be nonuniform, i.e., with increasing the density, 
the shape of nuclear matter region is changing from sphere 
(bec lattice) into cylinder, sl ab, cylindrical hole, and uni- 
form matter (inner fluid cor e) (jLorenz et al.fl993l ; IOvarnatsul 
1 19931 : ISumivoshi etaflll995h . This variation of nuclear struc- 
ture is known as the so-called "pasta structure" . The density 
that the cylinder structure appears, p p , depends on the nu- 
clear symmet ry energy expressed with the density symmetry 
coefficient L (|Ovamatsu fc Iidall2007n , which is suggested to 
be order p p ~ 10 13 g/cm 3 via the calculations of the ground 
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state of matter in the crust l|Lorenz et al.l 1 19931 ; lOvamatsul 
1 19931 ; ISumivoshi et~ahll 19951 ) . However, it might be quite dif- 
ficult to verify the existence of pasta structure by using the 
observation of neutron star properties such as mass and ra- 
dius, because the width of pasta phase is arou nd 10% of the 
crust , which is less than a few hundred meters ()Lorenz et al.l 
Il993h . In contrast, in this letter, we will calculate the tor- 
sional oscillations of neutron star with pasta phase and show 
that the frequencies of shear modes depend strongly on the 
presence of pasta phase and on p v . 

Previously, there are many ca lculations ab o ut th e 
shear modes ( e .g., [Leel (|2007h ; ISotani et all (|2007f ): 
ISteiner fc Watts! (2009D ), where they assume that nuclei 
form bcc lattice in the crust. With this assumption, the shear 
modulus of the crust is suggested as 

p = 0.1194 x m(Ze) 2 /a, (1) 

where is the ion number density, a — (3/47mi) 1//3 
is the average ion spac ing, and +Ze is the ion charge 
l|Strohmaver et all Il991). Roughly speaking, this relation 
can be expressed as th e power low with respect to the den- 
sity (|Sotani et alj|2007h . On the other hand, it is pointed out 
that the elastic properties in the pasta phase could be liquid 
cryst als rather than a crystalline solid (Pet hick fc Potekhinl 
1998). That is, unlike the case of bcc lattice, the shear modu- 
lus should be decreasing in the pasta phase as increasing the 
density. This picture is thought to be natural, because the 
structure of nuclear matter changes gradually as mentioned 
the above and at last the shear modulus becomes zero in the 
fluid core. In order to realize such a relation about the shear 
modulus, we adopt Eq. ((TJ in the crust except for the pasta 
phase, while in the pasta phase it is assumed that the shear 
modulus can be expressed as the cubic function with respect 
to the density, which satisfies that p should connect to Eq. 
|T} smoothly at p — p p and become zero smoothly at the 
boundary with the core, i.e., p = ci(p — p c ) 2 {p — C2), where 
ci and C2 are some constants determined by the boundary 
conditions (see Fig. [TJ. Of course, this simple relation of p 
in the pasta phase might be a kind of toy model, although 
expressed the rough behavior, because that should depend 
on the microscopic structure including th e matter composi- 
tion a nd/or the nuclear symmetry energy l|Ovamatsu fc Iidal 
I2007T ). But still, this relation is thought to be enough to ex- 
amine the dependence of shear oscillations on the existence 
of pasta phase as a first step. Anyway, as a result of the fall- 
off of shear modulus, one can expect that the shear velocity, 
v 3 = (p/p) 1 ^ 2 , will become smaller, and that the frequencies 
of shear modes will also decrease. In practice, although the 
understanding about the shear modulus of the liquid crystals 
is quite poor except for the suggestion that the shear mod- 
ulus becomes small in the pasta phase (|Pethick fc Potekhinl 
1998), it could be possible to obtain the information about 
the pasta phase and/or properties in liquid crystals via the 
observations of oscillations of neutron stars. Additionally, it 
should be noticed that the condition of p = at the crust- 
core interface corresponds to the limit of zero shear speed. 
Thus, the low-frequency modes might get trapped in the 
vicinity of crust-core interface and/or the large mode am- 
plitude might produce at the interface, which may in turn 
imply a large boundary-layer damping. 

To examine the shear modes, we prepare the static, 
spherically symmetric stellar models, which is the solution of 
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Figure 1. Shear modulus p as a function of relative radius r/R. 
Solid lines correspond to the shear modulus including pasta phase 
with different values of p p , while broken line is corresponding 
to that shown by Eq. (TTJ, where the stellar mass is fixed that 
M = 1.4M Q . 

the well-known Tolman-Oppenheimer-Volkoff (TOV) equa- 
tions described by a metric of the form 

ds 2 = -e^dt 2 + e 2A dr 2 + r 2 {d6 2 + sin 2 9d<t> 2 ), (2) 

where $ and A are the function of the Schwarzshild ra- 
dial coordinate r. Since the stiff equation of state (EOS) 
is favorable to the exp l anatio n of the observed QPO 
frequencies (ISotani et ah! 120071 '). we adopt the EOS L 
jPandharipande fc Smith! 1 19751 ) for the inner core. On the 
other hand, a modern EOS sugges ted in ( Douchin & Hacnsel 
l200lf ) is adopted for the crust (see lSotani et all l|2007h for the 
stellar properties), where the density at the boundary with 
core is set to be p c = 1.24 x 10 14 g/cm 3 l|Douchin fc Haensell 
120011 ). Then, the unknown parameter to determine the stel- 
lar model is only the density that the pasta phase appears, 
Pp- 

Since the pure torsional shear modes are incompressible 
and do not induce density variations in spherical stars, one 
can expect that no significant variation in the radiative part 
of the metric describing the gravitational field. Therefore, 
the frequencies of torsional shear modes are determined with 
satisfactory accuracy even when neglecting entirely the met- 
ric perturbations by setting <5<? M „ = 0, i.e., adopting the rel- 
ativistic Cowling approximation, and we adopt this approx- 
imation in this letter. The perturbation equation is derived 
from the linearized equation of motion, where one needs to 
provide the linearized shear stress tensor, 5T$ . In the same 
way as m ISchumaker fc Thornd (|l983l ). in this letter, 6T$ 
is assumed to be related to the linearized shear tensor, SS^, 
through ST^fJ — —2pSSfj,„. SS^v is determined from the re- 
lationship 8a nv = £ u 8S U v, where a uv d enotes the rate of 
shear tensor l|Schumaker fc Thorne|[l983l ). 

The torsional shear modes can be described with us- 
ing one perturbation variable, which is the angular displace- 
ment of the stellar matter, y(t,r). The non-zero component 
of perturbed matter quantities is oi-component of the per- 
turbed 4- velocity of fluid, <5w , which is expressed with y 
as 

Sv* = e -*d t ;y(i,r)-^<9eP,(cos0), (3) 

where dt and de denote the partial derivative with respect 
to t and 8, respectively, while Pi (cos 6) is the Legendre poly- 
nomical of order £. Assuming that the perturbed variable has 
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a harmonic time dependence, such that y(t, r) = e lult y(r) 
the perturbation equation reduces to 



y" + 



+ 



y 



- + - A' ) + — 



e+P 2 -2* (l + 2)(l-l) 
w e 



e A y = 0, (4) 



where e and p correspond to the energy density and pres- 
sure, respecti vely, and the prime denotes t he derivative with 
respect to r (jSchumaker fc Thorndll983t ). With appropri- 
ate boundary conditions, the problem to solve becomes the 
eigenvalue problem. We impose a zero traction condition at 
the boundary between the inner core and the crust (r = R c ), 
and the zero-torque condition at the stellar surface (r — R), 
which correspond to p,y' = at r = R c and y' _ at 
r = R (|Schumaker fc Thorndfl98^ ; ISotani et al.ll2007l ). So, 
if the shear modulus is exactly zero, the boundary condition 
at r = R c is automatically satisfied for arbitrary value of 
y' and the number of boundary conditions is not enough to 
solve the problem. However, since the shear modulus in the 
realistic case could be not absolutely zero but quite small 
value with the some fluctuation, in this letter we adopt that 
y' = even at r = R c . 

We examine the frequencies of torsional shear modes, 
as varying the value of p p in the range of p p = 10 13 — 10 14 
g/cm 3 , because that value is not certain but suggested to be 
order 10 13 g/cm 3 as mentioned before (|Lorenz et alJll993t 
lOvamatsul 1 19931 : ISumivoshi et al l Il995l ). Fig. [2] shows the 
fundamental frequencies of torsional shear modes with 1 = 2 
as a function of the stellar mass, where the broken line corre- 
sponds to the frequencies for the stellar model without pasta 
phase and the solid lines are corresponding to those with 
pasta phase with different values of p p . Obviously, one can 
observe that the frequencies of shear modes depend strongly 
on the existence of pasta phase. In fact, compared with the 
frequency for the stellar model without pasta phase, those 
with pasta phase are reduced to 12.0%, 34.9%, and 49.3% for 
p p = 1 X 10 14 , 4 x 10 13 , and 1 x 10 13 g/cm 3 , respectively. To 
compare with the observed frequencies, the lowest observed 
QPO frequency in the SGR 1806-20 (18Hz) is also plotted 
in this figure with the dot-dash line. From this figure, it is 
found that the smaller stellar mass is favored for smaller 
p p , if the frequency of 18 Hz is explained as usual with the 
fundamental I = 2 shear modes. Otherwise, with smaller p p 
than 6 x 10 13 g/cm 3 , one can explain the frequency of 18 Hz 
by using the fundamental shear mode with higher I. 

In Fig. [31 the frequencies of 1st overtones of shear modes 
with I — 2 are plotted as a function of the stellar mass. 
Unlike the fundamental modes, the frequencies with higher 
p p are almost same as that without pasta phase. This re- 
sult might be corresponding to the tendency in the Newto- 
nian limit, i.e., the frequency of fundamental modes could be 
roughly proportional to the shear speed, while those of over- 
tone depend on the ratio between curst thickness and stellar 
radius, which are almo st independent of the shear modulus 
IIHansen fc Cioffilll98ut ). Still, one can see the dependence 
of the frequency on the existence of pasta phase with lower 
pp. In fact, the frequencies with pasta phase are different 
from that without pasta phase in 0.4%, 46.0%, and 67.0% 
for p p = 1 x 10 14 , 4 x 10 13 , and 1 x 10 13 g/cm 3 , respectively. 
This dependence might be difficult to explain in analogy 
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Figure 2. Frequencies of fundamental torsional shear modes with 
I = 2 as a function of neutron star mass. The broken line cor- 
responds to the frequencies for the stellar model without pasta 
phase, while the solid lines correspond to those with pasta phase, 
where the labels denote the value of p p /(10 13 g/cm 3 ). Addition- 
ally, the dot-dash-line denotes the lowest observed frequency in 
the SGR 1806-20, which is 18 Hz. 
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Figure 3. Frequencies of first overtones of torsional shear modes 
with t = 2 as a function of neutron star mass, where the meaning 
of lines are same as in Fig. [2] The observed frequency in the SGR 
1806-20, which is 626.5 Hz, is also shown with the dot-dash-line . 



with the Newtonian limit, but one could be possible to ob- 
tain the additional information of the crust property via the 
observation of frequencies of overtones. Similar to Fig. [21 the 
observed frequency in the SGR 1806-20 (626.5Hz) is plotted 
in this figure with the dot-dash line, which is usually co nsid- 
ered to be caused by the 1st overtone of shear modes l|Pirol 
2005). But, with smaller p p , the frequency of 626.5 Hz might 
be corresponding to the 2nd or 3rd overtones. 

Since the both frequencies of fundamental and overtone 
shear modes depend strongly on the presence of pasta phase, 
whose effect has been neglected so far, one needs to consider 
this effect on the shear oscillations. Additionally, owing to 
this strong dependence, it could be possible to probe the 
properties of pasta phase via the observations of stellar os- 
cillations and stellar mass, although the constraint on the 
pasta phase is quite difficult via the other observations of 
neutron stars. In contrast to simple relation of p adopted in 
this letter, we will make an examination with more realistic 
shear model in the pasta phase in the future. 

At last, we will compare the calculated frequencies of 
shear modes with the observed QPO frequencies in giant 
flare. Especially, in this letter we focus on the QPO frequen- 
cies in the SGR 1806-20, i.e., 18, 26, 30, and 92.5 Hz in 
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Figure 4. Comparison of the frequencies of torsional shear modes 
(solid lines) with the observed QPO frequencies in SGR 1806-20 
(dot-dash-lines), where the adopted stellar model is with p p = 
1 X 10 14 g/cm 3 . 




In this letter, we consider the effect of nonuniform nu- 
clear structure, so-called "pasta", in the neutron star crust 
on the torsional shear modes. Based on the suggestion that 
the elastic properties in the pasta phase could be liquid 
crystals, the frequencies of shear modes are calculated with 
simple relation of shear modulus. As a result, due to the 
existence of pasta phase, one can observe the smaller fre- 
quencies than those expected without pasta phase. This 
result indicates not only the importance to take into ac- 
count the pasta phase, but also the possibility to probe the 
pasta structure via the observations of the stellar oscilla- 
tions, such as the QPO frequencies in giant flares. Further- 
more, we show the possibility to explain the observed QPO 
frequencies with using only crust shear modes with pasta 
phase, which is quite difficult with the traditional iden- 
tification without pasta phase. In order to find the most 
suitable stellar model with the observations, we need to 
adopt more realistic pasta structure and to examine sys- 
tematically with different EOSs, where the additional ef- 
fects , such as superfluidity inside the s tar (|Chamel fc Carter] 
l200fj ; iPassamonti fc Anderssoiill201ll ) as well as the stellar 
magnetic field, should be also taken into account. Still, we 
believe that the pasta structure in the neutron star curst 
could play an important role in the stellar oscillations. 

At the end, we are grateful to T. Tatsumi, K. Iida, K.I. 
Nakazato, N. Yasutake for their fruitful discussions, and also 
to reviewer for careful reading and giving the valuable com- 
ments. This work was supported by Grant-in- Aid for Scien- 
tific Research on Innovative Areas (23105711). 



Figure 5. Similar to Fig. [4] but with p v = 2 X 10 13 g/cm 3 . 

less than 100 Hz, because this phenomenon has the most 
obser ved frequencies among the detected giant flares in the 
past (|Watts fc Strohmaver|2006h . As mentioned in the intro- 
duction, some of the observed QPO frequencies are consid- 
ered as a result of the crust shear modes, while the others 
are associated with the magnetized fluid core. The reason 
is because of the difficulty to explain the all observed fre- 
quen cies with using only crust shear modes (|Sotani et al.l 
120071 ). That is, the fundamental £ = 2 mode, which is the 
possible lowest frequency, is considered to correspond to the 
observed frequency of 18 Hz and the fundamental £ = 3 
mode is corresponding to 26 or 30 Hz. But, the spacing of 
shear frequencies with different £ is larger than the spacing 
between the observed frequencies of 26 and 30 Hz. In prac- 
tice, according to such a traditional identification, we can 
explain the observed frequencies with the crust shear modes 
with pasta phase as shown in Fig. [4] i.e., 18, 30, 92.5 Hz can 
be identified as £ = 2, 3, and 10 fundamental modes within 
a few percent accuracy, where the expected stellar mass is 
M = 1.5M©. 

However, we find the possibility to explain the all ob- 
served frequencies with using only crust shear modes, if p p 
would be small. As shown in Fig. [5] the observed frequencies 
of 18, 26, 30, and 92.5 Hz can be identified as £ — 3, 4, 5, 
and 16 fundamental shear modes within a few percent accu- 
racy again, where the expected stellar mass is M = 1.5M©. 
This is important suggestion to explain the observed QPO 
frequencies in giant flares, which could become a directing 
post in the asteroseismology with neutron stars. 
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